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Abstract. We consider a free particle coupled with finite strength to a bath and 
investigate the evaluation of its specific heat. A harmonic oscillator bath of Drude 
type with cutoff frequency oJxj is employed to model an ohmic friction force with 
dissipation strength 7. Two scenarios for obtaining specific heat are presented. 
The first one uses the measurement of the kinetic energy of the free particle while 
the second one is based on the reduced partition function. Both descriptions yield 
results which are consistent with the Third Law of thermodynamics. Nevertheless, 
the two methods produce different results that disagree even in their leading 
quantum corrections at high temperatures. We also consider the regime where 
the cutoff frequency is smaller than the friction strength, i.e. ajjj < 7. There, we 
encounter puzzling results at low temperatures where the specific heat based on 
the thermodynamic prescription becomes negative. This anomaly is rooted in an 
ill-defined density of states of the damped free particle which assumes unphysical 
negative values when ')/u)y) > 1. 
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1. Introduction 

While some physical disciplines such as classical mechanics and electrodynamics 
underwent profound changes with the birth of quantum mechanics and relativity, 
thermodynamics proved impressively robust over the last century. The main reason is 
that the formulation of thermodynamics rests on few pillars only, such as entropy, 
temperature and the three Laws relating these state variables. The grandness 
of thermodynamics is that these concepts hold independently of the details of 
the corresponding total system dynamics. Nonetheless, the statistical mechanical 
foundation of thermodynamics strongly relies on the quantum mechanical properties 
of matter in particular what concerns the low temperature behaviour. 

Statistical mechanics gives rise to some subtle issues when going from a closed 
description of all degrees of freedom, including those of large environments, to a 
reduced description of an open system where all bath degrees of freedom are traced 
over. Generally no problems arise in the weak-coupling limit where the system-bath 
interaction tends to zero. This is the situation typically assumed explicitly, or at 
least implicitly, in the majority of textbooks. Even in this limit, however, pitfalls can 
arise in the quantum case, as recently elucidated in [1]. Generally, care must also be 
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taken in defining correctly the expression for (Gibbs) work in the First Law, as recent 
debates on the vaUdity of nonlinear fluctuation theorems have given evidence [2,3]. 

In this work we emphasize the finiteness of the coupling between the open 
quantum system and an environment of temperature ksT = 1/(3. The finite coupling 
will be modelled in terms of a generalized quantum Langevin equation (GLE) with a 
velocity-proportional memory friction, i.e. ohmic damping. An equivalent microscopic 
approach is provided by a bilinear coordinate-coordinate coupling between system 
degrees of freedom and environmental degrees of freedom, see e.g. in [1] and references 
given therein. 

The presence of an ohmic friction term leads to an irreversible motion with a 

unique stationary state. Already at this level, prominent differences with the classical 
situation emerge. As is well known, the classical canonical thermal equilibrium 
probability assumes the familiar Gibbs-Boltzmann expression oc exp{—0Hs), with 
Hs being the bare system Hamiltonian in absence of interaction. Most remarkably, 
this classical result holds true independently of friction strength. In striking contrast, 
the equilibrium density operator of an open quantum system becomes a function 
of the friction strength [1,4], thus exhibiting a dependence on the coupling to the 
environmental degrees of freedom. This feature can explicitly be inspected for an 
ohmic-like, damped harmonic quantum oscillator [1,5,6]: Only in the weak coupling 
limit does the canonical density operator reduce to the common Gibbs state, where 
Hs then represents the operator-valued system Hamiltonian. 

Given this observation in the preceding paragraph it should not come as too 
big a surprise that the evaluation of thermodynamic quantities for an open quantum 
system, such, for example, as its specific heat, is also plagued with subtleties. These 
difiiculties all originate from the assumed finite system-bath interaction. For the 
case of specific heat of a linear quantum oscillator of finite friction strength this was 
explicitly demonstrated in a recent work [7] . While in absence of damping zero entropy 
is approached exponentially fast, this approach is weakened with finite friction to a 
power law dependence in friction strength and temperature [7, 8] . The Third Law has 
also been validated for the quantum dissipative oscillator in the presence of velocity- 
coordinate and velocity-velocity couplings [9] as well as for a charged oscillator in a 
magnetic field [10]. The main finding is that the thermodynamic entropy of the open 
system vanishes according to a power law in temperature with the same exponent 
that characterizes the frequency dependence of the memory friction in the limit of 
vanishing frequency. 

It was remarked in [7], however, that the definition of the specific heat is 
ambiguous and may lead to pronouncedly diff'erent values. Therefore, the relation 
to experimental observations needs scrutiny, in particular in the study of nanosystems 
that behave quantum in nature but still are coupled with finite strength to an 
environment. Another surprising observation is that the familiar von- Neumann 
entropy for a quantum dissipative oscillator fails to approach zero for vanishing bath 
temperature [11]. 

In the following we elucidate in detail the complications that arise when evaluating 

thermodynamic quantities of open quantum systems whose dissipation strength is 
finite. A most suitable test bed is the case of a free quantum particle [4,7]. It is our 
working hypothesis that this simple system better be understood first before daring to 
embark on more complex physical situations. As is well known, a classical free particle 
does not obey the Third Law. Surprisingly, the coupling to a heat bath renders the 
system more quantum thereby helping to restore the Third Law [7] . 



The challenge of obtaining specific heat 



3 



The work is organized as follows. In the next section we couple a free quantum 
particle to a heat bath which exerts a finite dissipation on the free particle. We then 
focus on the evaluation of thermal equilibrium quantities, in particular the specific 
heat. In section [3] we introduce two definitions of the specific heat which a priori both 
seem physically well motivated. Their properties are investigated in sections |4] and 
m The emerging results cause worrisome "baffling" which we attempt to resolve by 
inspecting more closely the underlying density of states of a dissipative free particle. 
Some conclusions and consequences are given in our final section. 



2. Coupling a free particle to a heat bath 

2.1. Dissipation and corresponding quantum Langevin equation 

As mentioned in the introduction a free particle, or more generally, an ensemble of 
non-interacting, distinguishable particles fails to approach zero entropy at accessible 
temperatures. This classical behaviour is modified when the particle is able to 
exchange energy with a heat bath. We therefore focus on the simplest situation of a 
free quantum particle of mass M which is coupled to a heat bath made up of harmonic 
oscillators. The total Hamiltonian of particle plus bath thus reads (see figure [T]) 

H = Hs+Hb + Hsb (1) 

where 

describes the free particle with momentum operator p, and 

^ / 2 \ 

«B^i:(^+f.^f) (3) 

represents a set of harmonic oscillators constituting the heat bath which is bi- linearly 
coupled to the free particle via its position operator q by 

N 2 

i?sB = E^Hg., + ,2). (4) 

i=l 

Here, the coupling constants have been chosen without loss of generality in such a way 
that the ensemble of free particle and heat bath is translationally invariant so that 
the damped particle can still be considered as free. 

Elimination of the heat bath leads to an effective equation of motion for the 
position operator of the damped free particle, the quantum Langevin equation of a 
free particle, reading 

M^q + M j^ds^{t-s)^^q^m- (5) 
The coupling to the heat bath results in a damping kernel 
I ^ 

" M ^ '"^'^^^ cos(w,i) (6) 
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Figure 1. Setup of an open nanosystem with Hamiltonian Hs which is coupled 
with finite strength 7 through an interaction -ffsB to a bath described by a 
Hamiltonian Hg. The total system composed of system plus bath is weakly 
coupled to a super-bath which provides the working temperature T in thermal 
equilibrium. 



and an operator-valued noise term 

^{t) = -Mj{t - to)qito) 
N 

+ ^ [mi^tXiito) COS {uji{t~ to)) + uJiPiito) sin (uji(t- to))] , (7) 
1=1 

which depends on the initial conditions of free particle and heat bath. The sum of the 
second and third contribution of this expression, i.e. 

ij{t)^^{t)+M-f{t-to)qito), (8) 

only depends on the initial positions and momenta of the heat bath and is characterized 
by a vanishing first moment, 

m)B^o, (9) 

while the symmetrized correlation function obeys 



{r]{t)rj{s) + rj{s)r]{t))B = 3 2^ ^'^^ ^'^^ (^*(^ " \2k^) 



(10) 



2 

i=l 

Here, (. . .)b denotes an expectation value with respect to the equilibrium density 
matrix exp(— /JHb)/ Tr[exp(— /JiJe)] of the isolated bath. The temperature T is 
imposed on the heat bath via the weak coupling to a super-bath (see figure [ij . The 
commutator of the noise is non-vanishing, reading 

N 

[?]{t),7j{s)] = -ih^m^uj^ sin {uji{t - s)) , (11) 

and guarantees that the familiar commutator-relation between position q and 
momentum p is obeyed at all times, as required by quantum mechanics. 

In the following, we will restrict ourselves to the so-called Drude model of quantum 
dissipation where the damping kernel describes exponential memory on the time scale 
Wq^. For positive arguments t > 0, the damping kernel assumes the form: 

7(t) = 7WDexp(-WDi) (12) 

while for negative times we formally obtain ^{t < 0) = 7(|i|). This in turn assures 
that the symmetrized correlation function of the time- homogeneous noise correlation 
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in pUjl obeys a generalized Einstein relation [1]. The value 7 yields the damping 
strength according to 

di7(t)=7- (13) 





For later use we will also need the Laplace transform of the damping kernel, which is 
given by 

7(z) = / dtj{t) exp(-zt) = . (14) 

The total mass of the heat bath can be expressed via the appealing formula [4] 

y m, = M lim ^ . (15) 

i=l 

As a consequence of the fact that the Drude model behaves ohmic for low frequencies, 
i.e. 7(2:) goes to a positive constant 7 for z ^ 0, this heat bath behaves non-ballistic 
in the sense that it assumes an infinite mass. 

The Drude model is typically employed for regularization purposes. As an 
example we mention that in the case of strict ohmic dissipation, i.e. memoryless 
damping where j{t) = 2jS{t), the second moment {p^) of the momentum p with 
respect to the equilibrium density matrix of system plus bath exhibits an ultraviolet 
divergence. For the purposes of the regularization of such a quantity, it is assumed that 

is by far the shortest time scale in the problem apart possibly from the thermal 
time scale h(3. For a free particle subject to Drude damping one therefore typically 
assumes 7 <C wd- In this paper, we drop the requirement of a high-frequency cutoff 
which in turn enables us to study also structured, uncommon environments which 
may contain a low-frequency cutoff cjd, even smaller than the damping strength 7. 

The Drude model represents the simplest reservoir with memory in the sense that 
the memory effect can be described by a single additional degree of freedom. In fact, 
the deterministic equation of motion ([5]), i.e. where ^(t) is set to zero, is equivalent to 

q ^ V 

i) = z 

Z = —LOuZ — JLl!-£,V . (16) 

This system of differential equations contains a zero-frequency mode because system 
plus bath is translationally invariant. The remaining two eigenfrequencies are obtained 
from the second and third equation which describe a damped harmonic oscillator with 
damping strength wd and oscillator frequency (7Wd)^^^- For cou > 47, one finds 
exponentially damped motion which for very large cutoff frequencies contains the time 
scales 7 and lob as expected. However, for sufficiently small cutoff frequency, where 
WD < 47, the eigenvalues become complex and one observes a damped oscillation. 
This makes the Drude model more interesting than one might initially expect. 



2. 2. Thermodynamics of a dissipative free particle 

We note that a thermodynamic description of a free particle can be meaningful only 
if the particle is confined to stay within a box of finite size. The presence of the 
box leads to the quantization of the energy levels with the typical excitation energy 
AE = h^n'^ /2MLF' between the first excited state and the ground state for a one- 
dimensional box of length L with reflecting walls. This system in isolation (i.e. in the 
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absence of dissipation) approaches zero entropy only at extremely low temperatures 
when the thermal energy is comparable to or below the energy difference between the 
first excited and the ground state. 

We will restrict our investigations to temperatures which are high enough such 
that the discreteness of the spectrum of a single free particle can be ignored even in the 
weak coupling limit, i.e. we will always assume that PAE <C 1. For a Helium atom at 
a temperature of approximately 500 pK - the lowest temperature which can be reached 
today in the laboratory - the thermal energy and the excitation energy arc equal for a 
cavity of linear size of the order of 30 jjm. Even for such extremely low temperatures 
a handy box of linear size of, say, 1 cm would suffice to meet the required condition 
with l3AE « 10"-\ In order that any quantum effects can survive under thermal 
conditions rendering the discreteness of the particle spectrum practically invisible, 
the bath must provide relevant energy scales that are large compared to the thermal 
energy. In the case of a bath with a Drude cutoff, these relevant energies are given by 
the damping constant and the cutoff frequency. Hence, a regime of low temperatures 
with AE <C ksT <C h-y, hcori exists where quantum effects can be expected to become 
relevant. 



3. Free quantum Brownian motion: Two routes to calculate specific heat 

We are interested in the specific heat of a free damped particle. The volume in which 
the particle can move, will be assumed very large but fixed. In this case, the specific 
heat is obtained from the internal energy U by taking the derivative with respect to 
the temperature T, i.e. 

Usually, it is supposed that the coupling of the system to the heat bath defining 
the temperature can be treated in the limit of vanishing coupling strength. However, 
here we are interested in the case of finite coupling where the meaning of the system's 
internal energy is no longer obvious. In the following, we will study two different 
approaches. 

One possibility is to replace U by the energy E defined as expectation value of 
the system Hamiltonian Hg 

E = {Hs) , (18) 

where 

/u \ _ Trs+B[gsexp(-/3ij')] 

^^^^ - TrsMeM-PH)] ' ^ 
This leads to our first definition of a specific heat 

C^ = §- (20) 

This definition is based on the system Hamiltonian Hg and includes the interaction 
of the system with the bath only via the density matrix of the total system. 

Alternatively, one can start from the well-known and widely used expression for 
the partition function of the reduced system: It is defined in terms of the partition 
functions of the coupled system and of the uncoupled bath as [4, 5, 7, 8, 12-19] 

^ ^ Trs+B[exp(-/3g)] 

TrB[exp(-/3ifB)] ' ^ ' 
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where the total Hamihonian H consists of contributions from the system, the bath, 
and the couphng according to ([T]) . Employing the standard relation between partition 
function and internal energy 

U^-^HZ), (22) 

(j2ip implies that the internal energy is defined as 

U={H)- {Hb)b 

= E+ [{Hsb) + {Hb) - {Hb)b] ■ (23) 

While the internal energy U and the system energy E agree in the absence of a coupling 
between system and bath, this is no longer the case at finite coupling. Particularly 
disturbing is the observation that this difference is not solely given by the expectation 
value of the interaction Hamiltonian H^b , but still contains the difference of the bath 
energies caused by the interaction with the system. 

We can now proceed and define a second specific heat, reading 

The partition function (|2ip also allows us to define an entropy 



S = kT 



\n{Z) - (3^\nZ 



(25) 



This results from the thermodynamic relation 
dF 

S = -^ (26) 

with the free energy F — — (l//?)ln(Z). 

Alternatively, the second specific heat can as well be derived from the entropy 
(|25p: i.e. is equivalently obtained as 

C'=T§. (27) 

Our main concern is that these two routes of obtaining a specific heat generally 
yield different results. This will be elucidated further by studying the dissipation 
model of the Drude form when applied to the simple quantum dynamics of a free 
particle. 



4. Follow the route using energy 



We start our exploration of the specific heat by employing the definition (j20p based 
on the expectation value of the system energy 



(28) 



Here, M denotes the mass of the particle andp the momentum operator. For a general 
heat bath leading to an equation of motion of the form ([S]) and for thermal energies 
which are large compared to the energy scale T^-k^ j^Ml? one finds [4] 



1 



n=l 



(29) 
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Figure 2. Energy route: The specific heat is depicted as a function of 
dimensionless temperature for different ratios of ld-q/'-^. The approach towards the 
classical limit = 0.5 is at low temperatures inversely proportional to friction 
strength 7; it notably becomes enhanced for decreasing values of cutoff frequencies 
The specific heat itself is always positive and monotonically increases with 
increasing temperature T. 



with the Matsubara frequencies Un — 2Tm/h/3. In the high-temperature hmit, the sum 
vanishes and one recovers E = k^T /2, in agreement with the classical equipartition 
theorem. For strictly ohmic damping, 7(2) = 7, the sum diverges as mentioned above 
and one is obliged to introduce a cutoff in the spectral density of the heat bath. From 
(|29P we obtain by means of (|20|) the specific heat 



(30) 



where the prime denotes a derivative with respect to the argument. 

For strictly ohmic damping one has 7(2) = 7 and therefore ^'{z) — 0. It is 
tempting to drop the second term in the numerator of the sum and to evaluate the 
resulting sum which converges even for strict ohmic damping. However, this would 
yield a specific heat which diverges as the temperature approaches zero (see also the 
discussion after ([35]) ). The reason for this unphysical behaviour is that for constant 
7(2) the sum in (|29p does not converge and therefore derivatives should not be taken 
term-by-term. 

For the Drude model, one obtains with (fT4| for the specific heat C^: 



k-R 



X1X2 



where ^' {x) denotes the trigamma function and 



xi%l^'{xi)] - - , 



X1.2 



1± Wl- 



47 



(31) 



(32) 



47r \ V , 

The behaviour of this result is depicted in figure [5] for differently structured ohmic 
baths as characterized by various ratios of cutoff uj-q over damping strength 7. 

We next study some asymptotic limits more closely. For high temperatures T 
much larger than hj/k-g and huu/kB, the specific heat approaches the expected 
classical result with the leading quantum corrections reading 
_ 1 _ h^jivu 
"fe " 2 " 24(fcBT)2 



+ o(r-3). 



(33) 
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At low temperatures, i.e for k^T <C hj, hw-Q but still /cbT K^tt^ /2ML'^ the specific 
heat goes to zero linearly according to 



fee 3/17 15 V ^7 / V "^D/ 

This behaviour is in agreement with the Third Law: Finite quantum dissipation 
thus restores the validity of the Third Law of a free particle already at temperatures 
for which the discreteness of the particle spectrum is still irrelevant. The limiting 
behaviour is inversely proportional to damping strength. Put differently, strong 
dissipation diminishes the prefactor and thereby yields a turnover to the classical 
behaviour at higher temperatures only. The leading low-temperature behaviour of the 
specific heat only depends on the low-frequency behaviour of 7(2:). Notably, only the 
next-to-leading order also depends on the cutoff frequency wd- 

In the limit of an infinite cutoff frequency wd ~^ 00 the complete temperature 
dependence is obtained from (|3ip as 



fee \2TTkBTj ^ \2nkBTj 2nkBT 2' ^ ' 

The second term would have been missed in a naive calculation where one simply sets 
7'(z) = and subsequently evaluates the sum in (|30p . This term is important because 
it ensures that the specific heat does not diverge for T — > 0. 

5. Follow the route using the thermodynamic partition function 

In order to obtain the specific heat from the partition function (j2ip . one must first 
determine the latter quantity for the damped free particle. 

We again assume that the volume available to the particle is sufficiently large so 
that the discreteness of the energy levels can be neglected. In the limit of vanishing 
coupling to the heat bath, the partition fimction is given by 

7- /o \ 1/2 
L I 27rm ^ 



This result is a consequence of the density of states p{E) ^ E^^^^ of a free particle in 
one dimension. 

In the presence of a finite coupling 7 to the heat bath, the partition function is 
modified by the ratio of the two fluctuation determinants for the free particle in the 
absence and presence of the dissipative coupling [4,12,13,16]. We then obtain for this 
so defined partition function 

7-/0 \ 1/2 °° 

which will represent the starting point for all following calculations discussed in this 
section. We recall that the modifications due to the environmental coupling do not take 
into account the presence of confining walls. In using this expression for the partition 
function, we therefore assume a sufficiently large confining box so that boundary 
effects, which could show up in the partition function, can safely be neglected. 

The specific heat can be calculated from the partition function either by means 
of the internal energy (|24|) or the entropy (p7|) . Both routes lead to identical results. 
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In order to stay close to the reasoning of the previous section, we choose the first 
alternative. 

The internal energy of a damped free particle is obtained from (I37|) by means of 
reading 



U 



1 

2(3 



1 



n=l 



(38) 



This internal energy differs from the energy (j29p by the second term of the numerator 
appearing in the sum. 

For the Drude model, the sum can be evaluated and the internal energy is obtained 

as 



(39) 



where "^{x) is the digamma function and x\^2 are defined in (I32|) . This result should 
be contrasted with the energy E for which one finds from ([29|) for the Drude model 

1 



E = 



X\X2 



M N ('/'(a;i) - V'(a;2)) - „^ 



(40) 



For h(3u!-D/2Tr ;» 1, the internal energy U and the energy E thus differ by a constant 
energy contribution which does not play a role for the specific heat, but nevertheless 
indicates that the two quantities do not agree. This disagreement even persists in the 
case of strict ohmic damping. Note, however, that in this limit both quantities diverge 
logarithmically in the cutoff frequency wd, so that E and U will contain an infinite 
energy contribution independent of temperature. 

From ([59)1 it is straightforward to evaluate the specific heat by means of I 



(41) 



In the strictly ohmic limit lud oo, this result agrees with the corresponding 
expression ([35]) for . Some results obtained from (|4T|) are depicted in figure [3l 

For high temperatures but finite cutoff, i.e. if hP"f/2TT^ h(3uj-D/2Tr <^ 1, one obtains 

_ 1 _ n^juj-D 

fee 



o(r 



-3\ 



(42) 



! 12(fcBT)2 

While this reproduces the correct classical result, the leading correction differs by a 
factor of two from the high-temperature result of the specific heat ([55)1 derived from 
the system energy E. 

Very puzzling is the behaviour at low temperatures. In this limit the specific heat 
behaves like 



fee 



1 - 



1_ 



47r3 

Is" 



hrf 



1-3^- 



J_ 



O(T^). 



(43) 



Notably this result does again obey the Third Law in the same sense as does. 
Moreover, (|43|) agrees with ((34)) for wd ^ 7- The big surprise surfaces, however, for 
small cutoff frequencies. For lod < 7, the specific heat becomes negative at very 
low temperatures. This indicates a serious defect which must be hidden somewhere 
in the theory leading to (|43|). 
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Figure 3. Partition function route: The specific heat is depicted as a function 
of dimensionless temperature for different ratios of uiiy/'y. In the strict ohmic Hmit 
u)D — > this result agrees with , as depicted in figure [2] The behaviour as 
a function of finite cutoff is intriguing: While for ojd > 7i is monotonically 
growing with increasing temperature, as expected from thermodynamics, this fails 
to be the case when the friction strength surpasses the cutoff strength. Although 
still vanishing at absolute zero temperature we find at low temperatures formally 
negative- valued specific heats for this dissipative quantum system. 



The appearance of a negative specific heat is not restricted to the Drude model. 
This can be seen by expanding the general expression ([55]) for the internal energy in 
powers of temperature. By means of the Euler-MacLaurin summation formula one 
obtains 

2n-l 



U = Uo + Y^ 



B2n /27r 



^ (2n) 

11—1 ^ ^ 



/(2"-i)(oo)-/(2»-i)(0) 



(44) 



where i?2n are the Bernoulli numbers, /f*^' denotes the k-th derivative of / and the 
ground state energy is given by 



h r°° 

For the free damped particle, one finds by means of ([38]) 

7(3;) — X'j'{x) 



^^""^^ x + j{x) ■ 
From (I44p one readily obtains the specific heat 



(45) 



(46) 



fee 



B2 



27rfcBT 



-(2n-l)! 
The leading term of this expansion yields 
_ TT l + f (0) fceT ^ 

fee 



/(2"-i)(oo)-/(2"-i)(0)l . (47) 



(48) 



3 7(0) h 

in agreement with the specific heat ((43|) for the Drude model. This result implies 
that a negative specific heat will be found within the approach based on the partition 
function for every damping kernel with 7'(0) < — 1. 
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In contrast, for a damped harmonic oscillator with frequency loq, one obtains 

OJq + XJ(X) + 

resulting in the leading low-temperature behaviour of the specific heat 

Here, the specific heat remains positive independently of the damping kernel. 
6. Partition function and density of states 

The theory leading to the unexpected appearance of a negative specific heat, cf. 
(j43p . is based on an assumption and an approximation. The assumption is that 
the thermodynamic behaviour can be described by the partition function (j2ip . The 
approximation concerns the calculation of this partition function and limits the validity 
of the specific heat (H^ to sufficiently high temperatures for a fixed size of the box 
confining the particle. By choosing a sufficiently large box one though can approach 
such low temperatures for which the specific heat becomes negative but for which 
the approximation still is fully justified. 

It is therefore the form of the partition function as the ratio of the two partition 
functions of the total system and the bath which deserves further scrutiny. We will 
analyze the situation within the Drude model. Inserting the Laplace transform (|14p 
of the damping kernel into the general expression (|37p of the partition function, one 
finds 

L (2^mV'^ T{l + x^)T{l + X2) 

h\ [3 ) r(l + hPcuD/^n) ' ^ ' 

where r(z) is the gamma function and xi^2 are defined in ([5^ . 

Formally, the partition function can be related to a density of states p{E) of a 
damped system by means of a Laplace transform [17] 

/•OO 

Z{P)= dEp{E)exp{~f3E) . (52) 
Jo 

By its very definition as the number of states at energy E per unit energy the density 
of states p{E) must not assume negative values. This basic property of the density of 
states also restricts the admissible form of physically meaningful partition functions. 
In the remainder of this section we demonstrate that the unphysical behaviour of the 
specific heat as given by comes along with negative regions of the density of 
states. 

In order to discuss the partition function and the density of states it is useful 
to shift the origin of the energy scale to the ground state energy of the damped free 
particle 

Uo ^ ^In { + ^In { . (53) 

This expression is obtained from (j39p in the limit of zero temperature. Instead of the 
partition function Z wc thus consider Zexp(/3C/o)- 
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As the specific heat (j43p may become negative at low temperatures, we focus on 
the behaviour of the density of states at low energies. For ft/3wD 3> 1 one finds 

Zexp(/3C/o) = -^f^''^' 



(54) 



where 



is a characteristic length related to the Drude cutoff frequency ojd . The inverse Laplace 
transform of (|54p yields the low-energy behaviour of the density of states 

1/2 



piE) 



S 



E-U, 



.(56) 



Except for the delta function at the ground state energy Uq, the low-energy behaviour 
of the density of states is dominated by a constant which changes sign at 7 = lod . If 
the cutoff frequency lod is smaller than the damping constant 7, the density of states 
is no longer positive everywhere. This readily explains the surprising behaviour of the 
specific heat found in the previous section. In figure SI the energy dependence of the 
density of states is shown for wd/7 = 0.2, 1, and 5. For < 7, the density of states 
starts out with negative values as expected. At larger energies, peaks of the density 
of states are observed which can be related to the characteristic frequencies of the 
equation of motion (|16p which, for lod < 7 give rise to underdamped oscillations. For 
sufficiently large damping constant 7, there may exist several energy regions displaying 
a negative density of states. For 3> 7, one recovers the limit of vanishing damping 
where 

In the previous section, we have seen that the specific heat of a damped harmonic 
oscillator at low temperatures is always positive. An analysis along the lines presented 
above for the free damped particle indeed shows that the density of states, apart from 
a delta-function contribution, starts out at low temperatures with a positive value. 
However, the positivity condition of a density of states defined by (|52p is by no means 
obvious as we will discuss next. 

By its statistical mechanical definition the partition function is the sum over the 
Boltzmann factors counted with the multiplicity of the respective eigenenergies and 
therefore coincides with the Laplace transform of the density of states as a function of 
the energy. The corresponding Laplace variable is the inverse temperature /?. This in 
particular implies that the inverse Laplace transform of any proper partition function 
must not assume any negative values. 

Apparently, the positivity of the density of states is not automatically guaranteed 
by the above definition (PT|) as can be seen in the case of a free particle. Another 
example is provided by a system which is coupled to a rather trivial bath consisting of 
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Figure 4. Energy dependence of the density of states for uJo/'y = 0-2, 1, and 5. 
The soUd hne corresponds to the density of states of an undamped particle in one 
dimension. 



a single oscillator with angular frequency w. In this case one obtains for the partition 
function of the reduced system, according to the definition ((2T|) 

Trosciexp(/3i7osc)J n 

where En denotes the eigenenergies of the system coupled to the oscillator and the 
positive weights gn are the corresponding degeneracies. For a finite, spatially confined 
system, these eigenenergies constitute a discrete spectrum. The first factor of the 
right-hand side coincides with the partition function of the total system, Zs+b and 
the second term represents the inverse of the partition function of the bath consisting 
of a single oscillator. The inverse Laplace transform p{E) of Z becomes 

piE) = J2 9n5{E ~En + hj/2) - ^ gJiE - E^ - huj/2) (59) 

n n 

It consists of two sets of delta functions, one with positive weights at the energies 
En — huj/2 whereas the second set, which is located at the shifted energies En + huj/2, 
possesses negative weights. Therefore, neither does p{E) represent a physically 
meaningful density of states nor must its Laplace transform Z — J dE exp{— /3E)p{E) 
be interpreted as the partition function of the reduced system. 

For reservoirs which represent true thermal baths both the total system and the 
reservoir alone will have a continuous spectrum such that a more subtle compensation 
of the different contributions to the inverse Laplace transform may indeed lead to 
a positive density of states and a corresponding, physically meaningful partition 
function. 



7. Conclusions 



In this work we compared two definitions of specific heat for a free quantum particle 
which interacts at a finite strength with an oscillator heat bath possessing a Drude 
cutoff. The first definition was based on the assumption that the equilibrium 
expectation value of the Hamiltonian Hs of the system can be identified with the 
thermodynamic internal energy of the particle. The resulting specific heat approaches 
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the classical result for high temperatures and goes to zero for vanishing temperature 
in accordance with the Third Law. The second definition is based on the ratio of 
the partition functions of the total system and the heat bath which traditionally is 
identified as the partition function of an open quantum system [4,5,7,8,12-19]. 

The two definitions yield different results which even disagree in the high 
temperature regime in terms of their quantum corrections. If one associates (as done 
commonly) to each of the three partition functions of the total system, the system 
and the bath a free energy, the free energy of the open system is expressed as the 
difference of the free energies of the total system and the heat bath, in accordance 
with the "most remarkable formula" of Ford, Lewis and O'Connell [16]. We found 
that this definition may lead to unphysical results for the damped free particle, such 
as negative specific heat and negative regions of the density of states. A corresponding 
analysis of the respective quantities for a damped harmonic quantum oscillator though 
does not show any violations of their formal properties [7, 8] . This latter fact alone 
does not resolve the encountered baffling: The finding of a thermodynamic consistent 
result which is in no apparent violation with formal properties, such as a truly positive- 
valued density of states, does not necessarily imply that the so obtained values are 
physically meaningful. 

The definition of the specific heat via the expectation value of the system 
Hamiltonian Hg leads to a quantity that can be measured in principle. In contrast, no 
clear scheme seems to be available for a measurement of the specific heat that derives 
from the partition function of the reduced system. 

The conundrum of the physical meaning of a partition function of a reduced 
system that is coupled with finite strength to an environment thus still remains. 
We emphasize once more that the detected difficulties all "evaporate" when the 
dissipation can be described within the weak-coupling limit. However, the impressive 
experimental progress on tailored, well-defined nanosystems that behave quantum in 
nature but intrinsically are in interaction with omnipresent macroscopic environments 
requires to go beyond weak coupling. It thus calls for an urgent clarification of the 
thermodynamics of nanosystems and related issues. 
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